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Abstract. We report on Hartree-Fock-Bogoliubov (HFB) calculations of the ground-state rotational band 
of the heavy nucleus 254 No recently observed experimentally. The calculated quadrupole deformation is 
consistent with the experimental value of p = 0.27 and is almost constant over the whole band. We also 
reproduce fairly well the excitation spectra and moments of inertia of this isotope up to the maximal ex- 
perimentally observed state of spin 20. The rather high stability of this nucleus against fission is illustrated 
by the deformation energy curve providing very high fission barriers at zero spin within the HFB and HFB 
plus Lipkin-Nogami formalisms. The variation of these barriers with increased angular velocities is also 
studied. 

PACS. 21. 10. Re Collective levels - 21.60.Ev Collective models - 21.60.Jz Hartree-Fock and random- 
phase approximations 



After the Cohen, Plasil and Swiatecki seminal study of 
the stability of rotating liquid drops, very heavy nuclei with 
mass greater than 250 are generally expected to survive rota- 
tional excitation only up to rather low spins. Of course, large 
shell-correction energies added to the bulk liquid drop esti- 
mates may lead to sizable enhancements of the stability against 
fission in some cases. Observation of relatively high-spin states 
i in very heavy nuclei thus provides a very useful tool to assess 
any theoretical description of shell effects in this mass region, 
which is of primary importance for a good prediction of shell- 
stabilized superheavy elements. 

The experimental observation in the ground-state band of 
the 254 No nucleus of rotational states from spins 4 to 18/z (and 
tentatively 20/1) has been reported recently in a series of papers 
These experiments have been performed using 48 Ca 
beams from 130 to 219 MeV on 208 Pb targets. Two of them 
(including the more recent) were conducted at the Argonne 
National Laboratory with the Gammasphere 4jt Ge detector 
array, while the other made use of four Ge clover detectors 
and was performed at the University of Jyvaskyla. In all three 
experiments, the 254 No nuclei were implanted into a position- 
sensitive Si strip-detector placed at the focal plane. The y-rays 
coming from the 254 No were unambiguously identified from 
the coincidence with the oc-decay chain in the Si detector. As 
the transition to the ground-state have not been observed, the 
spin assignment relies on the use of a Harris parameterization 
for both /W and /' 2 ' moments of inertia. 

In ref. [Q], the authors proposed a new method to deduce the 
fission barrier height Bf based on a reconstruction of the entry 



distribution of the evaporation residues in spin and excitation 
energy. Since the fission process above the saddle-point energy 
should be much favored compared to the y-emission, they claim 
that for any given observed spin, the maximal excitation energy 
of the entry distribution should lie below the saddle-point en- 
ergy. For each spin the end-point of the excitation energy dis- 
tribution gives thus a lower bound for the fission barrier height 
which is estimated to be greater than 5 MeV around / = 12 h. 

In a recent paper involving some of the authors it has 
been shown using state-of-the-art semiclassical calculations, 
namely in the extended Thomas-Fermi (ETF) approach, that 
the rotating liquid drop model had to be refined by including a 
spin-dependence of its parameters which, as a result, provides a 
fission stability of heavy nuclei upon increasing the spin which 
is enhanced with respect to what is stated in ref. [jl]]. It is the aim 
of the present paper to go beyond this work, describing alto- 
gether bulk properties as well as shell and pairing effects within 
the fully self-consistent microscopic Hartree-Fock-Bogoliubov 
(HFB) formalism, in the particular case of 254 No. 

In our calculations, we have used the triaxial Hartree-Fock- 
Bogoliubov code developed by Laftchiev et al. [^,0]. This code 
is an extension of the Hartree-Fock Routhian code presented in 
refs. which uses Skyrme-type effective interactions and 
assumes parity and signature as symmetry operators for the 
one-body Hamiltonian. There, solutions breaking time-reversal 
symmetry are described through a decomposition of the single- 
particle wavefunctions on an axially symmetric harmonic os- 
cillator basis. The triaxial character inherent to the solutions of 
this formalism is taken care of by a decomposition of the vari- 
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ous densities as Fourier series in the azimuthal angle. This par- 
ticular choice provides shorter computation-times than usual 
triaxial codes as for instance those of refs. [ ^o[[ll[ ]. The pair- 
ing correlations are described a la HFB with particle-number 
projection tentatively taken into account within the approxi- 
mate Lipkin-Nogami (LN) scheme [ |i"2| , |l3| ] . The latter method 
which is widely used for instance in the context of fully self- 
consistent microscopic calculations (see e.g., ref. [ |l4| ] for time- 
reversal symmetric Hamiltonians, or ref. fll5[ l for Routhian cal- 
culations) results in the addition of a constraint on the second- 
order fluctuations of the particle-number. However it has been 
shown [|l6| in the A ~ 190 mass region from calculations using 
the Gogny force that this approximation could produce some- 
what inconsistent results. 

The HFB equations may be cast into a usual eigenvalue 
form in a doubled space where the 2x2 "Hamiltonian" is de- 
fined in terms of the Hartree-Fock Hamiltonian h and of the 
pairing potential A, whereas the eigensolutions are the usual 
U and V matrices (see for the notations, e.g, [|l7|]). Similarly 
to what was done in ref. jl5[], these equations are solved in 
two steps at each iteration. First we solve simple Hartree-Fock 
equations (that is the HFB equations with vanishing pairing 
field) to obtain the eigenstates of the Hartree-Fock-like Hamil- 
tonian. Then the original equations are solved in a truncated 
configuration space (through an energy cutoff) of the Hamilto- 
nian eigenstates previously determined. 

For the particle-hole channel, we have used the SkM* pa- 
rameterization of the Skyrme interaction fll8| ] since it was orig- 
inally fitted to provide a good description of the 240 Pu fission 
barrier and is thus rather well adapted to describe the stabil- 
ity against fission in very heavy nuclei. In the particle-particle 
channel, we have simply used the "modified seniority pairing 
force" previously described by Gall et al. in ref. [p5[]. At the 
present stage of our investigations of rotational properties in 
very heavy nuclei, we have not attempted to make a global 
fit over a variety of different nuclei to provide a pairing force 
which could be deemed to be somewhat universal in this mass 
region. We rather satisfied ourselves with a parameterization 
yielding the right moment of inertia at very low spins to as- 
sess the validity of our approach in reproducing the transi- 
tion energies (and moment of inertia at higher spins). The ma- 
trix elements G„ and G p for neutrons and protons are defined 
with the usual prescription G q = gq/(H +N q ), N q being the 
number of particles in the charge state q. The value of g„ (g p 
resp.) is 14.3 MeV (15.5 MeV resp.) in the pure HFB case and 
10.65 MeV (14. 1 MeV resp.) in the HFB+LN case with an en- 
ergy cutoff 6.2 MeV above the Fermi level in both cases. No 
smeared boundary conditions for this configuration space have 
been taken into account yet which may locally create some mi- 
nor technical deficiencies. 

The calculations have been performed using the usual de- 
formation-dependent truncation scheme equivalent to a 13 ma- 
jor shells spherical basis, corresponding to more than 450 har- 
monic oscillator states (for both signatures). The basis defor- 
mation parameters that have been used in our HFB calculations 
were obtained through a minimization of the energy calculated 
within an axially and time-reversal symmetric HF+BCS for- 
malism. At zero angular velocity non axially-symmetrical so- 
lutions have been forcefully searched out (due to the symme- 



try properties of the interaction, one would never explore as 
well known non-axial solutions starting from a first iteration 
axially-symmetrical ansatz). At finite spins, the solutions ex- 
plicitly break the axial symmetry. As a result however, the so- 
lutions are almost axially symmetric since the y angle never ex- 
ceeds 1° along the rotational band which allows us to only take 
into account the 5 lowest-order Fourier series in the densities 
decomposition (see ref. [|||). 

Recently, T. Duguet et al. Jl^ ] have presented the results 
of HFB+LN calculations which could be considered to be in 
essence rather similar to ours with two practical differences 
which are however not very substantial provided that the tech- 
nical work is performed adequately [(i) solutions of the HFB 
equations are computed on a spatial grid, (ii) a zero-range force 
is used in the particle-particle channel]. However they used 
for the particle-hole channel the SLy4 parameterization IpOlEl]] 
which is known to give in average (i.e., semiclassically) signifi- 
cantly too high fission barriers [E2[], This may be considered as 
a drawback to study fission stability properties. Another series 
of calculations for the same nucleus has been published quite 
recently [H making use of the D1S Gogny force parameter- 
ization [|24[ which is known to be well suited to study fission 
barriers. No corrections for the particle number symmetry vi- 
olation, e.g., of the Lipkin-Nogami type, have been taken into 
account. 

Along the rotational band of 254 No, the mass quadrupole 
moment Qo varies very slightly between 32.9 and 32.7 b in our 
calculations. This corresponds to a deformation parameter of 
P = 0.27 identical to the experimental value. It is to be noted 
that we are there in perfect agreement with the experimental re- 
sults [00,^1 as well as with the theoretical results of T. Duguet 
et al. [ [19[ yielding Q20 = 32.8 barn. The calculations of Egido 
and Robledo on the other hand yield slightly higher value of 
^20, resulting in a p value close to 0.29. 
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Fig. 1. Pairing energies (as defined in text) are plotted as functions 
of the angular velocity. Full and dashed lines (dash-dotted and dotted 
resp.) correspond to proton and neutron energies calculated within the 
HFB+LN (pure HFB resp.) formalism. 



We present in fig. |l] the evolution of the pairing energies, 
defined as a quantity proportional to the trace of the product of 
the abnormal density K with the pairing potential A, for protons 
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and neutrons with respect to the angular velocity Q.. It is well 
known from Nilsson-type calculations that for a deformation 
para meter of p ~ 0.3 the N = 152 neutron number is magic 
1 25]. This prediction is substantiated in our case by the rather 
low value of the neutron pairing energy compared to the pro- 
ton corresponding energy. As usual in HFB calculations, the 
pairing energy is seen to decrease with increasing angular ve- 
locity. In particular, the neutron pairing energy vanishes above 
h£l = 0.22 MeV. Within the same angular velocity range, the 
proton pairing energy is reduced by more than 40 percent. 

It is clear from fig. |], that around and above rotational fre- 
quencies of the order of h£l = 0.1 MeV the results of pure HFB 
calculations are rather dubious as far as the neutron pairing 
correlations are concerned. This remark is clearly applicable 
also to the calculations of Egido and Robledo. It is generally 
expected that the Lipkin-Nogami corrected HFB calculations 
should be better adapted to such a situation of weak pairing 
correlations. However, as already mentioned, the calculations 
of Peru et al. (See, e.g., JT^]) in the superdeformed A ~ 190 
region show that the quality of the results of such an approach 
is rather unpredictable. The application to such a Routhian ap- 
proach of pairing calculations conser ving ex plicitly the particle 
number, e.g., as those of Pillet et al. [ |26| , |27| ] would be of great 
interest there. In view of this expected inadequacy of pure HFB 
approach we have therefore limited the discussion below of cal- 
culated moments of inertia to HFB+LN results. 



are plotted in fig. ^| as functions of the angular velocity together 
with their experimental counterparts. 

The kinematic moment of inertia /W is, in fact, the expres- 
sion of the linear response of the system under the constraint 
of some collective variable. Nevertheless, we have actually as- 
sumed that it could be computed from the wavefunctions yield- 
ing the corrected energy (i.e., bypassing an explicit treatment of 
LN corrections specific to (j x ) in the spirit of ref. [|l4||). It may 
also be argued that due to the dependence of X2 upon the angu- 
lar velocity CO, equation (||) is not correct and dynamic moments 
of inertia should be evaluated through second derivatives of the 
energy with respect to the spin. However, we have checked that 
this co-dependence is very small and can be safely neglected (at 
least up to about CO = 0.1 MeV). 

It is seen on fig. ^ that we are reproducing rather well the 
experimental trend for both the /W and moments of in- 
ertia. The agreement of our results is better than the one ob- 
tained by Egido and Robledo for the only moment which they 
report, namely J^ K It is also better than the only moment re- 
ported by Duguet et al., namely in this case. The latter is 
quite remarkable because, as it can be seen on fig. ^[ our single 
particle spectra are somewhat different for the neutron states. 
Whereas these authors find important deformed gaps feature 
around N = 150 and N = 152 (see fig. 5 in [|l9|]), we get only 
a comparable N = 152 gap plus a N = 170 gap at low spins 
which is not present in the paper of Duguet et al. . 
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Fig. 2. Moments of inertia (see text) are plotted as functions of the 
angular velocity. Dynamic (kinematic resp.) moment of inertia calcu- 
lated within the HFB+LN formalism is displayed in full line (dashed 
resp.) while its experimental counterpart is represented by open circles 
(full resp.). 



The kinematic (/''') and dynamic (J^) moments of inertia 
are obtained from our calculations using the formulae 
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where j x is the total angular momentum operator and Q. is the 
angular velocity. The derivative involved in the calcula- 
tions have been taken over two consecutive even-spins. Both 
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Fig. 3. Single-particle Routhians in the No nuclei for neutrons (left 
panel) and protons (right panel) as a function of the angular velocity. 
The convention used for the (parity, signature) representation is the 
following: (+,+) in full lines, (+,-) in dashed lines, (-,+) in dash-dotted 
lines and (-,-) in dotted lines. Calculations are performed within the 
HFB+LN formalism. 



On fig. |] we display the fission barriers obtained at zero 
total angular momentum within the HFB and HFB+LN ap- 
proaches allowing triaxial deformation. The former yields a 
first fission barrier which is slightly lower than what is obtained 
for the latter. Comparing HFB+LN results now for two differ- 
ent total angular momenta / = Oh and / = 12 h we find that 
the first fission barrier height is slightly increased by 500 keV 
where one would expect it to be decreased. This phenomenon 
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is due to a higher single-particle level density at the top of the 
first barrier for / = 12/z as compared to what is obtained at zero- 
spin. At such low spins, shell structure effects are thus able to 
mask the anti-binding effect of the rotation leading to an extra 
fission instability With the basis size in use in our calculations 
we are not a priori able to provide reliable relative energies in 
the second fission barrier region, as noted years ago j 
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Fig. 4. Deformation energy curves are plotted as functions of the 
quadrupole moment Q20. Dashed line corresponds to HFB calculation 
at zero-spin, while full line (dotted resp.) corresponds to HFB+LN 
calculations at zero-spin (12h resp.). The energy reference is taken as 
the unconstrained zero-spin energy value in both formalisms. 



To conclude it appears that the SkM* parameterization of 
the Skyrme interaction (complemented by a seniority force in 
the p-p h-h channel) is able to provide a good reproduction of 
the spectroscopic data in the ground state rotational band of the 
254 No nucleus within the HFB plus Lipkin-Nogami approach. 
This result is all the more remarkable that similar calculations 
1 19 1 using a different effective force and yielding a rather dif- 
ferent shell structure in the relevant deformation range, lead to 
a similarly (may be slightly less) good reproduction of the ex- 
perimental data. Even though it is true that reproducing such 
collective rotation data provides a much wanted test of effec- 
tive forces to be used when predicting the stability properties of 
neighbouring superheavy nuclei, such a convergence of results 
from seemingly rather different effective forces should mitigate 
the hope to get so far here a completely stringent benchmark. 

The calculations of ref. [^3j have studied in detail the fis- 
sion properties of this nucleus, namely they underline its re- 
markable, and somewhat unexpected, fission stability which 
could partly be explained by a change of saturation properties 
of nuclear matter in presence of a centrifugal field as advocated 
in ref. |Q. Our results, limited in both spins and deformations 
as compared to those of ref. [^3| yield similar conclusions on 
this stability. A pending point common to all these calculations 
(those of refs. Jl9|p3|] as well as ours) is the pairing correlation 
treatment in weak pairing regions as encountered upon increas- 
ing the angular velocity. Surely, treatment of such correlations 
within an approach that explicitly conserves the particle num- 
ber, as the one proposed in refs.p(|^7||, would greatly improve 
the validity of the results. 



However whatever their mid stream character, the rather 
satisfactory reproduction of the difficult, and for that reason 
somewhat scarce, experimental results so far obtained by self 
consistent calculations as those of refs. [|l9j^3|] and the one pre- 
sented here, should constitute an incentive to experimentalists 
to provide more data so as to improve our knowledge of the ef- 
fective interaction to be used. This certainly constitutes a very 
timely endeavour. 
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